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Abstract 

We propose a CP-odd asymmetry in the supersymmetric process e + e~ — > XiXj ~^ 
XjT^f^ by means of the transverse polarization. We calculate the asymmetry 
and cross sections at a future linear collider in the 500 GeV c.m.s. energy range with 
longitudinal polarized beams and high luminosity. We work in the Minimal Super- 
symmetric Standard Model with complex parameters (jl, M\ and A T . The asymmetry 
can reach values up to 60%. We also estimate the sensitivity for measuring the r 
polarization necessary to probe the CP asymmetry. 



1 Introduction 



In supersymmetric (SUSY) extensions of the Standard Model (SM), some parame- 
ters can be complex. In the neutralino sector of the Minimal Supersymmetric Stan- 
dard Model (MSSM), these are the higgsino mass parameter \i and the gaugino mass 
parameter M±, while M<i can be chosen real by redefining the fields. In addition, in 
the scalar tau sector of the MSSM, also the trilinear scalar coupling parameter A T 
can be complex. The non-zero phases (p^, tpu^ an d <Pa t of these parameters give rise 
to CP-odd observables, which are not present if CP is maintained. Measurements of 
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such CP-odd observables will allow us to determine these phases, in particular also 
their signs. 



In this Letter we consider neutralino production 

e + e-^x-X?; i,J = l,...,4 (1) 

and the subsequent two-body decay of one neutralino 

X°^f±r T ; m = l,2, (2) 

for a fixed r-polarization. We would like to stress that without measuring the trans- 
verse polarization no sensitivity to the phase of A T , ifA T , can be obtained, because 
(2) is a two-body decay. When summing over the r~ polarization, we are sensitive 
only to CP violation in the production process [1,2]. The r~ polarization is given 
by [3] 

P " Tr(g) ' (3) 

with g being the hermitean spin density matrix of the r~ and <7j the Pauli matrices. 
We use a convention for P = (P!,P 2 ,P 3 ) where the component P3 is the longitu- 
dinal polarization and Pi is the transverse polarization in the plane formed by p e 
and p T . The component P 2 is the polarization perpendicular to p r and p e - and is 
proportional to the triple-product 

S r • (p r X p e -) , (4) 



where s r is the r~ spin 3-vector. Since under time reversal the triple-product changes 
sign, the transverse r~ polarization P 2 is a T-odd observable. Due to CPT invariance, 
P 2 is actually a CP-odd observable if absorbtive phases from final-state interactions 
are neglected. 

In this Letter we study the asymmetry 

A CP = l(P2-P2) , (5) 

which is CP-odd, even if absorbtive phases can not be neglected. In Eq. (5), P 
denotes the r~ polarization vector in the decay Xi ~^ ^m r_ an d P denotes the t + 
polarization vector in the decay xi ~^ ^ m T+ ■ I n Born approximation it follows from 
Eq. (5) that A C p = P 2 - 

In Section 2 we briefly review stau mixing in the MSSM and define the part of the 
interaction Lagrangian which is relevant for our analysis. In Section 3 we define the 
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t spin density matrix g and give the anlytical formulae. In Section 4 we discuss 
the qualitative properties of the asymmetry Aqp- We present numerical results for 
e + e~ — > Xi^iT m Section 5. We summarize and conclude in Section 6. 



2 Stau mixing and Lagrangian 



We give a short account of ft — tr mixing for complex fi = \/j,\e l(Pli , A T = \A T \e l<PAT 
and Mi = \Mi\e ltpM i . The masses and couplings of the r-sleptons follow from the 
hermitian 2x2 mass matrix which in the basis (tl,tr) reads [4,5] 



rr _ 



M? 

TLL 



e~ tipf \M? LR \ \ ( f L ^ 



e^\M? LR \ 



M? 

r RR 



(6) 



with 



M fLL =M L + (~0 + Sil12 COS2/9 ^ + 

M~ rr =M\- sin 2 9 W cos 2/5 m| + m 2 , 
Mf RL — (Mf LR )* — m T (A T — fi* tan (3), 
iff = arg[A r — ji* tan 0\, 



(7) 

(8) 
(9) 
(10) 



where m T is the mass of the r-lepton, ®w is the weak mixing angle, tan/3 = v-ijv\ 
with ^1(^2) being the vacuum expectation value of the Higgs field H^H®), and 
Mi, Mg,A T are the soft SUSY-breaking parameters of the fj system. The f mass 
eigenstates are (ri,r 2 ) = (tl,tr)7Z t with 



TV 



I e ilfif cos Of sin 6 f \ 
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(11) 
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The mass eigenvalues are 



< 2 = \ («, + M£J =F ^ML-M^ + ^M^p) . (13) 

The part of the interaction Lagrangian of the MSSM relevant to study the decay 
(2) reads (in our notation and conventions we follow closely [7,8]): 

C fTjc o = r k f(bl t P L + alP^x? + h.c. , % = 1, . . . , 4 , k = 1, 2 , (14) 

with 

4, = 9(K n y^ n , bl = g(Kl n yBl, (n = L, R) (15) 









( h T r \ 



















(16) 



and 

hli = W = YtN? 3 , 

fL = -^(ta,nQ w N a + N l2 ), 

r m = V2t a n<d w N* v (17) 

where Y T = m T /\/2m w cos /3, Pl,r = 1/2(1 =F 75) and g is the weak coupling con- 
stant. N is the 4x4 unitary neutralino mixing matrix, which diagonalizes the neutral 
gaugino-higgsino mass matrix Y a p, N* a Y a pNlp = m^oSik, in the basis (B, W 3 , H®, H®) 
[7]. The part of the Lagrangian for the neutralino production (1) can be found e.g. 
in [10,11]. 



3 Tau spin density matrix 

The unnormalized, hermitean, 2x2 spin density matrix of the r~ is defined by: 
g x " x 'k = |(|.M| 2 ) A ^dLips , (18) 

where M. and dLips are the amplitude squared and the Lorentz invariant phase 
space element, respectively, for the whole process of neutralino production (1) and 
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decay (2). The r helicities are denoted by A& and X' k . In the spin density matrix 
formalism (as used e.g. in [9,10]) the amplitude squared can be written as 

(\M\ 2 ) X ^ =2|A(x°)| 2 E M XtK (Pd)$- (19) 



It is composed of the unnormalized spin density matrices pp for the production (1) 
and po for the decay (2), the propagator A(%°) = l/\p Xi — ™> Xi + im x T Xi ], with 
p Xi ,m Xi , T Xi being the four-momenta, masses and widths of the decaying neutralino, 
respectively, pp and pp, carry the helicity indices \, A- of the neutralinos and/or the 
helicity indices A&, \' k of the r~. The factor 2 in Eq. (19) is due to the summation of 
the x°j helicities, whose decay is not considered. We introduce a set of spin basis vec- 
tors s" (a = 1, 2, 3) for the neutralino x? , which fulfill the orthonormality relations 
s Xi ' Xi = ~ anc ^ s Xi ' Pxi = 9- Then the spin density matrices can be expanded 
in terms of the Pauli matrices: 



(P 



p 



1 L ■}. 



P 5 



A; A'. 



Sa a 
P a Xi\ f •> 



(^)a* 



D XkX ' k 5 KXi + (£« ) x ^ aW 



(20) 



The analytical formulae of P and £ P can be found in [10]. Introducing also a set of 
spin basis vectors s b T for the r~, D XkX 'k and (E£,) AfcA fc can be expanded: 



D x ^=D5 XkK +D b a XkK , (21) 

(^l^-Syyi+S^ . (22) 

The expansion coefficient are given by 

D = Re{bl* al^m T m Xi + ^fc| 2 + |<J 2 )(p T • p Xi ) , (23) 

D b = lmMj 2 -\aU 2 )(PxrS b r)i (24) 

^D = l^xA\<f-\bl l \ 2 )(p T -sl), (25) 



^=Re(c;«D(^-4)(^-4) 

-(4 • ^)[^(lCl 2 + KJVr™*; + Re(b f m *a f mi )(p T ■ p Xi )} 
+Mbi:aUe^p T ^ Xiv sl/ T(7 , (26) 
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with e 0123 = 1. The last term in Eq. (26) contains the triple product (4). This term 
is proportional to Im(6^ a r m ^) and is therefore sensitive to the phases W and 
(Pmi- Inserting the density matrices of Eq. (20) into Eq. (19) yields 

(\M\ 2 ) X ^ = 4|A(x°)| 2 [(PD + ^ a D )5 XkK + {PD b + X a P V$)a b XkX , )] . (27) 



4 Transverse tau polarization and CP asymmetry 

From Eqs. (3) and (27) we obtain for the transverse r~ polarization 
n J |A(x°)| 2 dLips 



/ |A(x°)| 2 PP> dLips 



(28) 



which follows because in the numerator we have used / |A(x°)| 2 PD 2 dLips = 
and in the denominator we have used / |A(x°)| 2 SpEJ, dLips = 0. As can be seen 
from Eq. (28), P 2 is proportional to the spin correlation term S^ 2 , Eq. (26), which 
contains the CP-sensitive part lm(b'^*a^ ni )e fMUpCT p T ^p^. v s'i^s 2 a . In order to study the 
dependence of P2 on the parameters, we expand 



lm(bl*a f u ) = /[cos 2 9 f Y T lm(f T u N l3 ) + sin 2 f Y T V2 tan Q w lm(N a N i3 ) 

+ sm 2 9rCos 2 9 f (Y 2 lm(N l3 N i3 e llpf ) + V2 tan G^mi^A^e-^))] , (29) 

using Eqs. (15)-(17) for m — 1. If CP violation is solely due to (fiA T 7^ (mod it), 
we find from (29) that P2 oc sin 26 f sin ip f . We note that the dependence of tpf on 
ipA T is weak if \A T \ <C \fi\ tan/3, see Eq. (10). Thus, we expect that P 2 increases with 
increasing \A T \. 

Details concerning phase space and kinematics necessary for the calculation of P 2 
from Eq. (28) can be found in [2]. The r~ spin vectors s b are chosen by: 

1 f n s 2 x s 3 \ 2 / p T x p e - \ 3 1 / E T \ 

S r= '| 1 ' S r= M 1 > S r = Pr , 1 fPr • 30 

V S 2 XS 3 / V Pr X Pe- / m T \ p r I 



In order to measure P 2 and the CP asymmetry .Acp, Eq. (5), the r~ from the decay 
(2) and the r + from the subsequent f + decay, f+ — > x? r+ 5 have to be distinguished. 
This can be accomplished by measuring the energies of the r's and making use of 
their different energy distributions [2]. 

A potentially large background may be due to stau production e + e _ — > Tj + f ~ — > 
T+T ~XiXv However, these reactions would generally lead to "two-sided events", 
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whereas the events from e + e~ — > XiXi ~^ r+r XiXi are "one-sided events". More- 
over, the background reaction is CP-even and will not give rise to a CP asymmetry, 
because the staus are scalars with a two-body decay. 



5 Numerical results 

We present numerical results for e + e _ — > X1X2 an d the subsequent decay of the 
neutralino into the lightest stau x\ ~^ f° r a linear collider (LC) with y/s = 500 
GeV and longitudinal polarized beams with (P e -, P e +) = (0.8, —0.6) or (P e -, P e +) = 
(—0.8,0.6). This choice favors right or left selectron exchange in the neutralino 
production process, respectively [10]. We study the dependence of the asymmetry 
.Acp and the production cross sections a = a p (e + e~ — > X1X2) x BR(x2 ~^ ^i T ~) 011 
the parameters (p^, </?Mi, \Mi\, (Pa t , and tan/3. 

For the calculations we assume \Mi\ = 5/3 tan 2 Q W M 2 , m T = and use in Eqs. (7) 
and (8) the renormalization group equations (RGEs) for the selectron masses [6], 
M\ = m 2 + 0.79M| and M| = mg + 0.23M|, taking m = 100 GeV. The size of \A T \ 
is restricted due to the tree-level vacuum stability conditions [12]. The restrictions 
on the masses of the SUSY particles are m-± > 104 GeV, > 100 GeV and > 
m^o. For the calculation of BR(%2 — > t~\T ~) we concentrate on the parameter domain 
where two-body decays are allowed and neglect three-body decays. We consider the 
two-body decays 

X°2 - r m r, £ RtL £, x° 2 Z, jfiW* xlHl £ = e,fi, m, n = 1, 2, (31) 

with H® being the lightest neutral Higgs boson. The Higgs mass parameter is 
chosen as itla = 1000 GeV, which means that explicit CP violation is not important 
for the lightest Higgs state [13]. Furthermore, the neutralino decays into charginos 
and the charged Higgs bosons X2~/> X.t^ T y as we ^ as decays into the heavy neutral 
Higgs bosons X2~/> Xi -^2,3 > are ruled out in this scenario. 

In Fig. 1 we show the contour lines for Aqp hi the <y9^ T -|A T | plane. As can be seen 
Aqp is proportional to sin26'f simp^, which is expected from Eq. (29). „4cp increases 
with increasing \A T \ 3> |//| tan/3, which also follows from Eq. (29). Furthermore, in 
the parameter region shown the cross section a varies between 20 fb and 30 fb. 

In Fig. 2 we show the dependence of Aqp on tan (3 and (pu 1 ■ Large values up to ±20% 
are obtained for tan/3 « 5. Note that these values are obtained for p Ml ~ ±0.87r 
rather than for maximal pu x ~ ±0.57r. This is due to the complex interplay of the 
spin correlation terms in Eq. (28). In the region shown in Fig. 2, the cross section o 
varies between 10 fb and 30 fb. 

Figs. 3a and 3b show, for <pa t = 0.57T and <p Ml — ip^ — 0, the \n\-M 2 dependence of 
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Fig. 1. Contour lines of „4cp in Eq. (5) 
for M 2 = 200 GeV, \n\ = 250 GeV, 
tan P = 5, = = an d 

(P e _,P e+ ) = (0.8,-0.6). 
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Fig. 2. Contour lines of „4cp in Eq. (5) 
for A T = 1 TeV, M 2 = 300 GeV, 
= 250 GeV, ipA T = (p^ = and 
(P e -,P e+ ) = (0.8,-0.6). 



the cross section o and the asymmetry .Acp, respectively. The asymmetry reaches 

values up to —15% due to the large value of \A T \ = 1 TeV and the choice of the beam 

polarization (P e -,P e +) = (—0.8,0.6). This choice also enhances the cross section, 

which reaches values of more than 100 fb. The gray shaded area excludes chargino 

masses m^± < 104 GeV. In the blank area either the sum of the masses of the 
1 

produced neutralinos exceeds \/s = 500 GeV or the two-body decay ~^ ^i r ^ is 
not open. 

For ip Ml = 0.57T and ip^ = <pa t = we show in Figs. 4a,b the contour lines of a and 
^4cp in the |/z|-M 2 plane, respectively. It is remarkable that in a large region the 
asymmetry is larger than -10% and reaches values up to -40% while also the cross 
section is large. Unpolarized beams would reduce the largest values of a by a factor 
3, whereas Aqv would only be marginally reduced. 

For \fi\ = 300 GeV and M 2 = 400 GeV, we show in Figs. 5a,b contour lines of a and 
Acp, respectively, in the (p ll -(fM 1 plane. As can be seen the asymmetry Aqp is very 
sensitive to variations of the phases tpM 1 and ip^. Even for small phases, Acp can be 
sizable. Small values of the phases, especially of </? M , are suggested by constraints on 
electron and neutron electric dipole moments (EDMs) [15] for a typical SUSY scale 
of the order of a few 100 GeV (for a review see, e.g., [16]). 

The polarization of the r is analyzed through its decay distributions. The sensitivi- 
ties for measuring the polarization of the r lepton for the various decay modes are 
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Fig. 3. Contour lines of a and -4cp in the |//|-M 2 plane for ipA T = 0.57T, (fM 1 = ^ = 0, 
A T = 1 TeV, tan/3 = 5 and (P e -,P e +) = (—0.8,0.6). The blank area outside the area 
of the contour lines is kinematically forbidden since here either y/s < + m^ 2 or 
mfj + m r > m£ 2 . The gray area is excluded since m~± < 104 GeV. 
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Fig. 4. Contour lines of a and -4cp in the |/^|-M 2 plane for ifM 1 = 0.5tt, (fA T = <£n = 0, 
A T = 250 GeV, tan/3 = 5 and (P e -,P e +) = (-0.8,0.6). The blank area outside the 
area of the contour lines is kinematically forbidden since here either y/s < + m^ 2 or 
mfj + m T > m^ 2 . The gray area is excluded since m~± < 104 GeV. 
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(j(e + e — > Xi T i T ) i n fb 





Fig. 5. Contour lines of <r and »4cp i n the f^-fMi plane for M2 = 400 GeV, = 300 
GeV, tan/? = 5, (p Ar = 0, A T = 250 GeV and (P e _,P c +) = (-0.8,0.6). 

quoted in [17]. The numbers quoted are for an ideal detector and for longitudinal r 
polarization and it is expected that the sensitivities for transversely polarized r lep- 
tons are somewhat smaller. Combing informations of all r decay modes a sensitivity 
of S — 0.35 [18] has been obtained. Following [17], the relative statistical error of 
P 2 (and of P 2 analogously) can be calculated as 5P 2 = AP 2 /|P 2 | — cr s / (-SI-P2I y/~N), 
for a s standard deviations, where N = aC is the number of events with integrated 
luminosity C and cross section a = a p (e + e~ — > X1X2) x BR(x 2 ~~ *■ ^i T ~)- Then for 
Acp, Eq. (5), it follows that A^4 C p = AP 2 /v / 2- We show in Fig. 6 the contour lines 
of the sensitivity S = \/2/(\Acp\y/N) which is needed to measure .Acp at 95% CL 
(cr s = 2) for C =500 fb -1 , for the parameters (fA T = 0.2n, ifu-i = V 9 ^ = 0- m Fig. 7 
we show the contour lines of the sensitivity s for the parameters ^>m x = 0.27T and 
¥v = { Pa t = 0. It is interesting to note that in a large region in the \fi\-M 2 plane in 
Figs. 6 and 7 we obtain a sensitivity S < 0.35, which means that the asymmetries 
can be measured at 95% CL. 



6 Summary and conclusion 



We have proposed and analyzed the CP odd asymmetry Acp in Eq. (5) in neu- 
tralino production e + e~ — > x^X^j an d the subsequent two-body decay of one neu- 
tralino Xi ~^ ^ r=F - The asymmetry is due to the transverse r T polarization, which 
is non-vanishing if CP-violating phases of the the trilinear scalar coupling param- 
eter A T and/or the gaugino and higgsino mass parameters M±, \i are present. The 
asymmetry occurs already at tree level and is due to spin effects in the neutralino 
production and decay process. In a numerical study for e + e~ — > X1X2 an d neu- 
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Fig. 6. Contour lines of S for ipA T = 0.2tt, 
= ¥7, = 0, A T = 1 TeV, 
tan/3 = 5 and (P e -,P e +) = (-0.8,0.6). 
The blank area outside the area of the 
contour lines is kinematically forbidden 
since here either y/s < + m^ 2 or 

mfj +m T > m^ 2 . The gray area is excluded 
since m~± < 104 GeV. 
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Fig. 7. Contour lines of S for 99^1 = 0.27T, 
9?a t = ^ = 0, A T = 250 GeV, 
tan/3 = 5 and (P e -,P e +) = (-0.8,0.6). 
The blank area outside the area of the 
contour lines is kinematically forbidden 
since here either y/s < + m^ 2 or 

mfj +m T > . The gray area is excluded 
since m-± < 104 GeV. 



tralino decay % 2 ~" * r 1 ± r T we have shown that the asymmetry can be as large as 
60%. It can be sizeable even for small phases of \x and M 1; which is suggested by the 
experimental limits on EDMs. Depending on the MSSM scenario, the asymmetry 
should be accessible in future electron-positron linear collider experiments in the 500 
GeV range. Longitudinally polarized electron and positron beams can considerably 
enhance both the asymmetry and the production cross section. 
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